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We report on our observations of the parsec-scale radio jet structures of five blazars that have been detected by 
ground-based TeV gamma-ray telescopes. These five blazars all belong to the class of High-frequency peaked 
BL Lac objects (HBLs), which are the most common blazar type detected at the TeV energy range. Because 
of their relative faintness in the radio, these HBLs are not well represented in other radio blazar surveys. Our 
observations consist of five epochs of Very Long Baseline Array (VLBA) imaging from 2006 to 2009, of each of the 
five blazars 1ES 1101-232, Markarian 180, 1ES 1218+304, PG 1553+113, and H 2356-309, at frequencies from 
5 to 22 GHz. Fundamental jet properties, including the apparent jet speeds, that can be measured from these 
multi-epoch series of VLBA images are presented and compared with other gamma-ray blazars. Confirming 
prior work, we find that the TeV HBLs have significantly slower apparent jet speeds than radio selected blazars. 
Together with other radio properties of the HBL class, this suggests modest Lorentz factors in their parsec-scale 
radio jets. This in turn suggests some form of Lorentz factor gradient in these jets, since they are likely to 
have high Lorentz factors in their TeV-emitting regions. The study presented here approximately doubles the 
number of TeV HBLs with multi-epoch parsec-scale kinematic measurements. 
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1. INTRODUCTION 

The new generation of ground-based TeV gamma- 
ray telescopes has now detected about 50 AGN as 
of this writing (tevcat.uchicago.edu). Most of these 
belong to the class of high-frequency peaked BL Lac 
objects (HBLs), with synchrotron and inverse Comp- 
ton spectra peaking in the X-rays and high-energy 
gamma-rays, respectively. These blazars display re- 
markable variability in their gamma-ray emission on 
time scales of minutes (e.g., [!(), suggesting extremely 
high Lorentz factors of about 50 in the TeV emitting 
region of their relativistic jets (e.g., Q). 

The only way to image these relativistic jets directly 
on the parsec-scale is in the radio with VLBI. Most 
HBLs are relatively faint in the radio, so they are 
not included in other flux-limited VLBI monitoring 
projects. We are currently studying the parsec-scale 
structures of the TeV HBLs with multi-epoch VLBI 
imaging using the National Radio Astronomy Obser- 
vatory's Very Long Baseline Array (VLBA). Our re- 
sults for six of these TeV blazars have already been 
reported [|[, here we present multi-epoch VLBA ob- 
servations of an additional five TeV blazars discovered 
by the TeV telescopes during 2006. 

Studying the parsec-scale structure of these radio- 
weaker high-frequency peaked BL Lac objects with 
VLBI is fundamentally important to blazar physics, as 
they are likely to represent a physically distinct AGN 
class from the flat-spectrum radio quasars and the low 
and intermediate- frequency peaked BL Lac objects. It 
has been suggested that the former group (the HBLs) 
represent non-evolving low-excitation radio galaxies 
(LERGs) with less efficient or absent accretion disks 
and weak jets, while the latter group (the FSRQs, 
LBLs, and IBLs) represent evolving high-excitation 



radio galaxies (HERGs) with standard accretion disks 
and powerful jets Q. It has also been suggested 
that unified schemes applied to the weak-jet sources 
are consistent with velocity gradients in the radiating 
plasma in the jets of these weak-jet sources and 
such velocity gradients may be directly observable in 
VLBI imaging in the form of transverse structures in 
the images. 



2. OBSERVATIONS 

We observed the five TeV blazars 1ES 1101-232, 
Markarian 180, 1ES 1218+304, PG 1553+113, and H 
2356-309 at five epochs each from 2006 to 2009 with 
the VLBA under observation code BP146, at frequen- 
cies from 5 to 22 GHz, for a total of 25 new VLBA 
images produced for this paper. With the VLBA, ob- 
serving at higher frequencies provides superior angular 
resolution but worse sensitivity. For each source ob- 
served for this paper, the observation frequencies were 
chosen to obtain the best possible angular resolution 
while still significantly detecting the parsec-scale jet. 

Most observations recorded three hours of data 
on source, at a data rate of 256 Mbps. Polariza- 
tion observations were made of Markarian 180 and 
PG 1553+113 (the two brightest sources) at the first 
epoch. The data were calibrated in AIPS, and im- 
aged in Difmap. All sources were detected with direct 
fringe fitting at all epochs, with peak flux densities 
ranging from 14 to 148 mjy beam -1 . Phase refer- 
encing was not needed for delay and rate corrections; 
although phase-referenced observations were made of 
the three sources 1ES 1101-232, 1ES 1218+304, and 
H 2356—309 at the first epoch in order to determine 
more accurate milliarcsecond scale positions of these 
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Figure 1: VLBA images of 1ES 1101-232, Markarian 180, 1ES 1218+304, PG 1553+113, and H 2356-309 at 8 or 
22 GHz from 2009 Jul 2. The axes are labeled in milliarcseconds (mas). The lowest contour in each image is three 
times the rms noise level in that image, and each successive contour is a factor of two higher. Peak flux densities are 
23, 60, 29, 91, and 22 mjy beam -1 for Figures a to e, respectively. The rms noise levels are 0.05, 0.22, 0.04, 0.14, and 
0.05 mjy beam -1 for Figures a to e, respectively. Positions of Gaussian component centers are indicated by diamonds. 
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sources than were previously available. 

An image of each of these five blazars, from the 
2009 Jul 2 epoch (the final epoch), is shown in Fig- 
ure 1. After imaging, we fit circular or elliptical Gaus- 
sian components to the calibrated visibilities using the 
modelfit task in Difmap. Note that fitting to the visi- 
bilities allows sub-beam resolution to be obtained, so 
that components are clearly separated from the core in 
the model fitting even when they appear blended with 
the core in the images. The complete set of 25 images 
and the corresponding models will be published else- 
where. After model fitting, linear fits were made to 
the separations from the core of the Gaussian model 
component centers versus time, in order to measure 
the apparent jet speeds. For this fitting, error bars 
were assigned to the model component positions us- 
ing the method described by @. 



Table I Apparent Component Speeds 



Proper Motion Apparent Speed 
Source Comp. (mas yr" 1 ) (multiples of c) 
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3. THE INDIVIDUAL SOURCES 
3.1. 1ES 1101-232 

Observations of this blazar at a redshift of z=0.186 
were made at frequencies of 5 and 8 GHz. Its structure 
consists of a 25 mJy core, and a short jet extension 
to the northeast that is represented as a single 5 mJy 
component in our modeling. 

3.2. Markarian 180 

Markarian 180 (1133+704, z=0.045) is among the 
two brightest of these five blazars in total radio flux 
density, and we observed it at a frequency of 22 GHz. 
At this frequency, a jet extends several mas to the 
east from a 40 mJy core. This jet is modeled by four 
Gaussian components at most epochs. VLBI polar- 
ization observations were done at the first epoch, and 
they showed that the VLBI core is 4% polarized, with 
the electric vector position angle aligned with the jet 
position angle. 

3.3. 1ES 1218+304 

Observations of this z=0.182 blazar were made at 
frequencies of 5 and 8 GHz. The images show a jet 
with a large apparent opening angle extending about 
10 mas to the east from a 25 mJy core. The inner 
3 mas of this jet are modeled as two stationary Gaus- 
sian components in our observations. 

3.4. PG 1553+113 

The redshift of this source is somewhat uncertain. 
The redshift has been constrained to be in the range 
from 0.43 to 0.58 Q, and we have used the mean of 



this range, or z=0.5, for results requiring a redshift 
measurement. This source is relatively bright, and 
our VLBA observations were made at a frequency of 
22 GHz. This source is quite compact at 22 GHz, with 
a ~ 100 mJy core and a short jet to the northeast that 
can be modeled as a single ~ 25 mJy component about 
0.2 mas from the core. A single epoch of polarization 
observations at the first epoch showed the core to be 
3% polarized, with the electric vector position angle 
aligned with the jet position angle. 

3.5. H 2356-309 

Observations of the z=0.165 source H 2356—309 
were made at 5 and 8 GHz. H 2356—309 is both the 
faintest and the southernmost of the five blazars ob- 
served in this series, making the imaging more diffi- 
cult. The parsec-scale structure of this source consists 
of a ^20 mJy core and a jet extending several mas to 
the southeast that is modeled by two components. 



4. RESULTS 

For these conference proceedings, we present re- 
sults on the apparent jet speed measurements in these 
five blazars. Other results from this series of VLBA 
observations, including transverse structure measure- 
ments, polarization images, brightness temperatures, 
and flux variability measurements will be presented in 
a more complete paper elsewhere. 

All together, we measured apparent speeds for 10 
different jet components in these five blazars, ap- 
proximately doubling the number of TeV HBLs for 
which multi-epoch parsec-scale kinematics informa- 
tion is available. These new apparent speed measure- 
ments are tabulated in Table 1. Apparent speed values 
were calculated assuming cosmological parameters of 
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Figure 2: Distribution of fastest apparent jet speeds in the 11 TeV HBLs studied by us to date (top), compared with 
the fastest apparent jet speeds measured in the radio-selected blazar sample from the Radio Reference Frame Image 
Database pj (bottom). All apparent speeds are measured from first order linear fits to the separations from the core of 
the Gaussian model component centers versus time. For the TeV HBLs, six of the speeds are from Table 5 in while 
the remaining five come from this paper. Formally negative speeds, most of which are statistically consistent with no 
motion, are included in the leftmost bin. 
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H Q = 71 km s" 1 Mpc" 1 , Q m = 0.27, and fL A = 0.73. 
Note that there is a large range in the proper mo- 
tion error associated with the different components. 
In general, the relatively bright compact components 
detected at all epochs in the high-frequency observa- 
tions have more precise proper motion measurements 
than the fainter more diffuse components detected in 
the lower-frequency observations. 

The apparent speeds presented in Table 1 are all 
low, with values ranging to only slightly over lc, and 
the measurements actually show no detected motions 
at all that are above 2er significance. This is in marked 
contrast to radio and GeV-selected blazar samples, 
which have typical apparent jet speeds of several c 
(e.g., 0, H). The distributions of the fastest mea- 
sured apparent jets speeds in each jet for the eleven 
TeV HBLs that we have observed to date and for a 
radio-selected blazar sample are compared in Figure 2. 
The distributions differ at >99.99% confidence accord- 
ing to a Kolmogorov-Smirnov (KS) test. This differ- 
ence between the HBLs and the more powerful blazars 
was previously noted for the first six TeV blazars stud- 
ied with the VLB A [3J, and now these additional five 
TeV HBLs confirm and statistically strengthen that 
earlier result. 



5. DISCUSSION 



factor gradients have also been proposed to explain 
the Lorentz factor differences between the radio and 
TeV-emitting regions in these sources, see H and ref- 
erences therein for a more complete discussion of mod- 
els in the literature. 

Of the additional newly discovered TeV HBLs that 
have been detected by the TeV telescopes within the 
past couple of years, a number have quite low flux den- 
sities in the radio, below that of the blazars that were 
imaged for this paper. However, the sensitivity up- 
grades recently completed at the VLBA should allow 
the parsec-scale structures of all of the extragalactic 
TeV detections to be imaged in the near future, and 
definitively compared with the more powerful sources 
detected by Fermi that are already well-studied in the 
radio. 
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When considered together with several other lines 
of evidence discussed in detail in 3 and [HI, this lack 
of superluminal apparent speeds in the TeV HBLs sug- 
gests that the Lorentz factor is relatively low in the 
parsec-scale radio jets of these sources (r ~2-3) , which 
contrasts with the high Lorentz factors that are typ- 
ically required to explain their TeV gamma-ray vari- 
ability and spectral energy distribution. 

These modest radio-jet Lorentz factors imply that 
the relativistic jets of the TeV HBLs may be struc- 
tured, with a Lorentz factor gradient either along the 
jet (so that the jet decelerates by the parsec scale, 
e.g., [Lll) or transverse to the jet (a fast spine sur- 
rounded by a slower layer or sheath, e.g., 12]), or 
both. The existence of limb brightening in the VLBI 
images of the nearest TeV blazars (Markarian 421 and 
501, 3, O, [3, [H) supplies additional direct evi- 
dence that transverse structures exist in these jets. A 
number of other models in addition to simple Lorentz 
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